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Abstract 

A measurement of the inclusive deep inelastic neutral current e^p scattering cross section 
is reported in the region of four-momentum transfer squared, 12GeV^ < < 150 GeV^, 
and Bjorken x, 2 ■ 10"^ < x < 0.1. The results are based on data collected by the HI Col- 
laboration at the ep collider HERA at positron and proton beam energies of E^, = 27.6 GeV 
and Ep = 920 GeV, respectively. The data are combined with previously published data, 
taken at Ep = 820 GeV. The accuracy of the combined measurement is typically in the 
range of 1.3 - 2%. A QCD analysis at next-to-leading order is performed to determine the 
parton distributions in the proton based on HI data. 
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1 Introduction 



The electron-proton collider HERA extends the kinematic range of deep inelastic lepton-nucleon 
scattering (DIS), determined by the four-momentum transfer squared, Q^, and Bjorken by two 
orders of magnitude towards high and small x, compared with fixed target experiments. This 
allows proton structure to be thoroughly investigated as is vital for the understanding of strong 
interactions, described by Quantum Chromodynamics (QCD). 

This paper presents the most accurate cross section data to date for the inclusive neutral 
current process e^p e^X, measured in the kinematic region 12GeV^ < < 150GeV^ 
and 2 • 10""^ < x < 0.1. The data newly presented here were taken in the year 2000 with 
positrons of energy Eg = 27.6 GeV and protons of energy Ep = 920 GeV, corresponding to a 
centre of mass energy ^/s = 319 GeV. The luminosity amounts to 22 pb ^ The measurements 
are combined with similar data taken in 1996/97 at Ep = 820 GeV [1]. The present paper 
refers to methods and detailed explanations given in a recent publication [2] of lower data, 
0.2 GeV^ < < 12GeV^. The results can be compared with data obtained by the ZEUS 
Collaboration [3,4]. 

The double differential neutral current DIS cross section in its reduced form and neglecting 
contributions from Z boson exchange, is given by 

with or denoting the fine structure constant and f{y) = y^/[l + (I - y)^]. The inelasticity y 
is related to Q^, x and the centre-of-mass energy squared, s = AE^Ep, hy y = I sx. A 
first measurement of the longitudinal structure function at low x was recently performed by 
HI [5]. The measurement presented here is restricted to the region of inelasticity y < 0.6 where 
the contribution of Fi is small and thus it focuses on the structure function F2. The data are 
used to determine the derivative (5^2/(9 In Q^)^ which provides a sensitive test of the evolution 
dynamics of partons, and the derivative {dXnFjIdXn x)qi, which quantifies the rise of F2(x, Q^) 
at fixed towards low x [6] . 

The new data cover the region of deep inelastic scattering, from a few GeV^ to about 
150 GeV^, with unprecedented accuracy at low x. A QCD analysis at next-to-leading order 
(NLO) is performed to obtain a new set of parton distribution functions (PDFs) from the inclu- 
sive DIS cross section measurements of the HI experiment alone. The QCD analysis is based 
on the results presented here, the low data [2], and the neutral and charged current (NC and 
CC) data sets at high [7-9]. 



2 Measurement Technique 

The analysis techniques are similar to those applied at lower [2], more details of the present 
analysis can also be found in [10]. 
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2.1 Detector 



The HI detector [11,12] was built and upgraded for the accurate measurement of inelastic ep in- 
teractions at HERA. The detector components most relevant to this measurement are the central 
tracker, the backward' lead- scintillator calorimeter (SpaCal) [13] and the liquid argon calorime- 
ter (LAr) [14]. The central tracker consists of the central jet drift chamber, two complementary 
z drift chambers, the central inner (CIP) [15,16] and outer proportional chambers, and the cen- 
tral silicon tracker (CST) [17]. The drift chambers and the CST are used for the measurement 
of tracks from the hadronic final state. The momenta of the tracks are determined and the event 
vertex is reconstructed. The polar angle of the scattered positron is determined by the planar 
backward drift chamber (BDC) and the vertex position. Complementary tracking information 
is obtained from the backward silicon tracker (BST) [18]. The SpaCal contains electromag- 
netic and hadronic sections. Its energy resolution for electromagnetic energy depositions is 
6EIE « 0.07/ V-E/GeV © 0.01. It also provides a trigger based on the scattered positron en- 
ergy. The LAr allows the hadronic final state to be reconstructed. Its energy resolution was 
determined to be SE/E « 0.50/ ^/EJGeV © 0.02 with pion test beam data [19]. Two electro- 
magnetic crystal calorimeters, a photon tagger and an electron tagger, located close to the beam 
pipe at z = -103.1 m and z = -33 m, respectively, are used to monitor the luminosity via the 
measurement of the Bethe-Heitler process ep yep. 

2.2 Online event selection 

The online trigger conditions used in this analysis are based on an energy deposition in the 
electromagnetic section of the SpaCal. Three trigger conditions with different energy and ra- 
dius thresholds are used, which largely overlap in phase space. The dominant source of trigger 
inefficiency are the veto conditions against beam related background. For these a global ineffi- 
ciency of (0.5 + 0.3)% is determined and corrected for. The inefficiency of the online software 
filter is determined to be 0.2%, which is applied as a global correction with a systematic uncer- 
tainty of half that size. The residual trigger inefficiencies from other sources are smaller than 
0.1%. 



2.3 Kinematics 

The DIS event kinematics are reconstructed from scattered positron and hadronic final state 
information. The positron energy E'^ and scattering angle 9e are used. For the hadronic final 
state a sum over the particles' energies Ej and longitudinal momenta P^j is formed, S/, = YiiiEj- 
P^i). The total difference between energy and longitudinal momentum E - P^ h obtained by 
adding to Z/, the positron contribution, E - P^ = l,^ + E'^{1 - cos 6^). Based on these variables, 
measurements of and y are obtained using the electron and the S methods as explained 
in [2]. In order to optimise the measurement accuracy, the electron method is used at larger 
y > OA, while the S method is used at lower y. For this analysis, positrons with scattering angles 

'The backward direction is determined by the outgoing positron beam direction. HI uses a coordinate system 
with the positive z axis given by the outgoing proton beam direction and the nominal interaction point at z = 0. 
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between 6, « 158° and Oe « 173°, energies E'^ > 11 GeV and Ql > lOGeV^ are included. The 
cross section measurement is performed in bins of x and chosen similarly to the previous 
measurement [1], with small modifications due to the different centre of mass energy. For each 
bin and reconstruction method the purity P = N^scgen/Nrec and the stability S = A/^rec,gen/A^gen are 
calculated. Here A^rec (A^gen) is the total number of reconstructed (generated) Monte Carlo events 
in the bin and A^recgen is the number of events which are both generated and reconstructed in the 
same bin. The values for purity and stability exceed 40% in all analysis bins and are typically 
well above 50% for the chosen reconstruction method [10]. 

2.4 Positron and hadronic final state reconstruction 

The reconstruction of the scattered positron is based on the measurement of a deposition of 
energy, termed a cluster, with a limited transverse size characteristic of an electromagnetic 
shower. The energy of the cluster E'^ is obtained by summing over all cells of the cluster in 
the electromagnetic section of the SpaCal. Its transverse size is characterised by Riag, which 
is obtained from the positions of all SpaCal cells belonging to a cluster using a logarithmic 
energy weighting [20]. A cut Riog < 4 cm is applied. For additional background suppression, 
the energy deposition in the hadronic section of the SpaCal behind the electromagnetic cluster. 
Eh, is required to be less than 15% of E'^. The positron candidate cluster is further required to 
be associated to a track in the BDC, formed by at least 4 hits from the 8 layers. This ensures an 
accurate measurement of the polar angle 6^ in combination with the z position of the interaction 
vertex, Zvtx, determined with the central track detectors. 

The reconstruction of the hadronic final state uses information from the central tracker and 
the calorimeters LAr and SpaCal [2] . The determination of S/, is affected by the presence of 
noise in the calorimeters. The resulting bias is particularly strong for small = S/,/2£'g. 
Contributions of noise from the SpaCal and the LAr are suppressed as described in [2]. 

2.5 IMonte Carlo event simulations 

Monte Carlo (MC) simulations are used to correct for detector acceptance and resolution effects, 
and for the background subtraction. The DIS signal events are generated using the DJAN- 
GOH [21] event generator, while the PHOJET [22] program is used for the photoproduction 
background. Elastic QED Compton events are generated using the COMPTON event gener- 
ator [23]. The cross section measurement is corrected for QED radiation up to order a using 
HERACLES [24]. The radiative corrections are cross checked with HECTOR [25]. An agree- 
ment to better than 0.3% is found in the kinematic range of this measurement. 

All generated events are passed through the full GEANT [26] based simulation of the HI 
apparatus and are reconstructed using the same program chain as for the data. For consistency, 
the calibrations of the SpaCal and the LAr, as well as the alignment, are performed for the 
reconstructed MC events in the same way as for the data. The calorimeter noise is determined 
using events from random triggers and is overlaid on the simulated events. The simulated 
events are reweighted to match the cross section derived from the QCD fit which is described 
in section 6. 
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Description 



Requirement 



Kinematic Range 

Scattered positron energy 

SpaCal cluster radius 

Energy in hadronic SpaCal section 

BDC validation 

Radial cluster position 

Vertex z position 

Transverse momentum balance 

Longitudinal momentum balance 

QED Compton Rejection 



Ql > lOGeV^ 
E'^>n GeV 
i?iog < 4 cm 
EJE'^ < 0.15 



> 4 linked hits, BDC-SpaCal radial match < 2.5 cm 



''spac < 73 cm 
Izvtxl < 35 cm 



P^/Pip > 0.3 
E-P-> 35 GeV 
Topological veto 



Table 1 : Event selection criteria. 



3 Data Analysis 



The data for the present measurement were recorded in the year 2000 and correspond to an 
integrated luminosity of 22 pb"' . In the following a description of the analysis is given. Further 
information can be found in [2,10]. 

3.1 Event selection 

An overview of the selection criteria is given in table 1 . Detection of the scattered positron is 
required within the backward calorimeter SpaCal, as discussed in section 2.3 and 2.4. In order 
to ensure that the cluster is well contained in the SpaCal, the extrapolation of the associated 
BDC track segment to the SpaCal plane is required to be within a distance from the beam line 
of ?"spac < 73 cm. The event vertex is reconstructed either by the central drift chambers and the 
CST or by using the CIP and the positron cluster position. Its z-position is required to be within 
+35 cm of the centre of the detector. 

Events for which the hadronic final state is poorly reconstructed are rejected by demanding 
that the total measured hadronic transverse momentum Pj be at least 30% of the positron trans- 
verse momentum P^. This efficiently removes migrations from very low y. Events with high 
energy initial state photon radiation are excluded by requiring E - > 35 GeV. The QED 
Compton process is suppressed using a topological cut against events with two back-to-back 
electromagnetic clusters reconstructed in the SpaCal. 

3.2 Efficiency determination 

The efficiencies of the positron identification requirements (transverse cluster size, hadronic en- 
ergy fraction, BDC validation) exceed 99% in most of the phase space. They are well described 
by the simulation. The only significant local inefficiency, of about 5%, is observed for the BDC 
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validation requirement at a radial distance from the beam line of tbdc ~ 25 cm where the geom- 
etry of the BDC drift cells changes. A detailed map of this inefficiency as a function of tbdc and 
the azimuthal angle of the scattered positron is obtained using positron candidates validated by 
the BST. The simulation is adjusted accordingly and an additional 0.5% systematic uncertainty 
is added to the cross section measurement error, uncorrected from bin-to-bin. 

The efficiency of the vertex reconstruction is determined using events independently recon- 
structed by the BST. This efficiency is determined to be close to 100% for all but a few bins at 

< 20 GeV^ and y < 0.03, where it drops to about 75 - 95%. It is described by the simulation 
to an accuracy of 0.3% for y > 0.01. For lower y the description is accurate to about 1%, which 
is accounted for by an additional uncorrelated uncertainty in the corresponding bins. 

3.3 Alignment and calibration 

The alignment of the HI detector starts from the internal adjustment of the central tracker and 
proceeds with the backward detectors, BDC, SpaCal and BST. The alignment of the BDC is 
performed using the tracks of positron candidates that are reconstructed in the central tracker 
with high accuracy. The SpaCal position is adjusted based on the positron tracks measured in 
the BDC. Finally the BST is aligned using events with a well reconstructed central vertex and a 
positron track measured in the BDC. The resulting agreement of the polar angle measurements 
is better than 0.2mrad, which defines the associated systematic uncertainty. The measured and 
simulated scattered positron polar angle distributions are shown in figure la). The data are well 
described by the simulation. 

The calibration of the electromagnetic scale of the SpaCal corrects for differences in the 
gain factors of individual SpaCal cells, for local non-uniformities at the sub-cell level and for 
global non-linearity in the energy response. The calibration is based mainly on the positron can- 
didates at low y. For these the kinematics can be reliably reconstructed using the double angle 
reconstruction method as described in [2], which employs only the polar angle information of 
the hadronic final state and the scattered positron, in addition to the positron beam energy. The 
non-linearity of the energy response is determined and corrected for using a sample of yy 
events. The energy scale is then checked using elastically produced 7/i/r particles decaying to 
e^e and QED Compton events, ep epy, where the scattered positron and photon are both 
reconstructed in the SpaCal. The relative data-to-simulation scale uncertainty is derived from 
two contributions: a global 0.2% arising from the double angle calibration and a part resulting 
from the studies, which is 1% at £ = 2GeV linearly decreasing to zero at E = 27.6 GeV. 
The measured and simulated scattered positron energy distributions are in very good agreement, 
as shown in figure lb). 

The calibration of the calorimeters employed for the hadronic final state energy measure- 
ment is based on kinematic constraints relating the scattered positron to the hadronic final state. 
For the calibration of the LAr calorimeter, the conservation of the total transverse momentum 
Pt is used. The hadronic SpaCal calibration utilises the conservation of E - P,. 

The transverse momentum balance between the scattered positron and the hadronic final 
state is studied as a function of various variables, such as Pj, the polar angle of the hadronic 
final state, and y^. Figure Ic) shows the overall P'r^/Pf^. distribution with a vertical line at 0.3 
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indicating the analysis cut value. For Pj/Pj values larger than this cut, the data distribution 
lies inside a band given by varying the LAr hadronic energy scale in the simulation by 2% 
for y > 10~^. At the lowest y ^ 0.005 considered in the measurement, the hadronic final 
state is produced at small polar angles, and a significant part escapes detection. The systematic 
uncertainty on the hadronic energy scale is increased from 2% at y = 10"^ linearly in logy to 
10% at y = 10'^. Furthermore, for bins with y < 10"^, an additional uncorrelated cross section 
uncertainty of 2% is estimated by varying the Pl^/Pj cut between 0.25 and 0.35. 

A topological finder is used to identify and subtract LAr noise. The fraction of hadronic 
energy attributed to noise is described by the simulation to within 15%, which is taken as a 
systematic uncertainty. 

For large values of y, the contribution of the SpaCal to S/, becomes larger than the combined 
contribution of the LAr calorimeter and the tracks, and thus the total E - P,, expected to be 
2Ee, provides the calibration for the hadronic final state measurement in the SpaCal. A study 
at y ^ 0.5 shows that the hadronic energy measurement in the SpaCal is described by the 
simulation to 0.3 GeV. Figure Id) shows the E - P. distribution for the data and the simulation. 
The simulation reproduces the data within the combined calibration uncertainties. 

3.4 Background 

A small source of background for this analysis arises from photoproduction events, in which 
the scattered positron escapes detection in the backward beam pipe while a particle from the 
hadronic final state mimics the positron. For a fraction of photoproduction events the scattered 
positron is detected by the electron tagger of the luminosity system. The photoproduction MC 
(PHOJET) is normalised globally based on a comparison with tagged events applying all selec- 
tion criteria apart from the E - Py cut. The systematic uncertainty on the background normal- 
isation is taken to be 15%, based on the studies described in [2]. The background contribution 
at the highest y ~ 0.5 considered here typically amounts to 5%, hence yielding an uncertainty 
of less than 1% on the cross section. Potential background from non-ep interactions is studied 
using non-colliding bunches and found to be negligible. 

3.5 Summary of systematic uncertainties 

The systematic uncertainties are classified into two groups, bin-to-bin correlated and uncorre- 
lated errors. They are summarised in table 2. The electromagnetic and hadronic energy scales, 
the positron scattering angle, LAr noise, background subtraction and normalisation are all con- 
sidered to be correlated sources of uncertainty. The uncorrelated errors arise from various effi- 
ciencies and the radiative corrections. For most of the analysis phase space none of the sources 
of systematic uncertainty dominates the result. For very low y < 0.01, many uncertainties in- 
crease strongly, as the reconstruction efficiency drops rapidly and the modelling of the hadronic 
final state suffers from increased losses in the forward direction. 

The uncertainty on the global normalisation of the measurement is determined mostly from 
the luminosity measurement, which is accurate to within 1.1%. The additional corrections 
discussed in section 2.2 increase the overall normalisation uncertainty to 1.2%. 
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Correlated systematic errors 



Source 



Uncertainty 



E'^ scale 
E'^ linearity 
Polar angle 6e 
LAr hadronic scale 

LAr noise contribution io E - 

SpaCal hadronic scale 

yp background normalisation 

Luminosity and other global uncertainties 



0.2% 

1% at 2 GeV to 0% at 27.6 GeV 

0.2 mrad 

2% for i/ > 0.01 
increasing linearly in log y to 10% at «/ = 0.001 

15% 

0.3 GeV 

15% 

1.2% 



Uncorrelated systematic errors 



Source 



Uncertainty 



BDC efficiency 
Vertex efficiency 
Radiative corrections 

Additional uncertainty for bins with y <0.0\ 



0.5% 
0.3% 
0.3% 
2.0% 



Table 2: Summary of the systematic uncertainties. For the correlated error sources, the uncer- 
tainties are given in terms on the uncertainty on the corresponding source. For the uncorrelated 
error sources, the uncertainties are quoted directly in terms of the effect on the measured cross 
section. 
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3.6 Data - Monte Carlo Comparison 



The quality of the MC description of the measurement can be verified by comparing experimen- 
tal and simulated distributions in the kinematic range of the measurement, corresponding ap- 
proximately to j/s > 0.005 and ye < 0.6. The DIS MC cross section prediction is reweighted to 
the QCD fit discussed in section 6. The distributions are normalised to the luminosity. Figure 2 
a)-d) shows the x and distributions of the selected events for the two kinematic reconstruc- 
tion methods, electron and S. A good overall agreement is obtained in the description of the 
data by the simulation. 



4 DIS Cross Section Results 

The measurement of the inclusive double differential cross section for deep inelastic positron- 
proton scattering, e^p e^X with the Ep = 920 GeV data, is reported in the tables 6 and 
7. The tables show the statistical, uncorrelated and the various correlated error contributions. 
The region of small Bjorken x, 2 ■ 10~^ < x < 0.1, and four-momentum transfer squared, 
12GeV2 <Q^ < 150 GeV^, is covered. 

The stability of the cross section measurement is tested with a set of dedicated cross checks. 
This is done by splitting the data into two approximately equal subsamples and comparing 
the results obtained using the same reconstruction method on each. For example, the data 
are compared as measured with the upper and the lower half of the SpaCal, for negative and 
positive Zvtx positions, and dividing the sample into early and late data taking periods. These 
tests are sensitive to local efficiency problems, energy miscalibrations, or the stability of the 
luminosity measurement. A particularly interesting test is the comparison of the cross section 
measurements performed with the electron and E methods, which have different sensitivities to 
systematic error sources. The results shown in figure 3 demonstrate very good agreement taking 
the correlated uncertainties into account. For the final result, the method with the smaller total 
uncertainty is chosen in each bin, which results in a transition near y = OA. 

This analysis improves the uncertainties by up to a factor of two with respect to previously 
published results on the inclusive DIS cross section in this kinematic range by HI [1] and by 
ZEUS [3,4]. 



5 Combination of 820 and 920 GeV Data 

The present measurement is based on data taken with 920 GeV proton beam energy. A similar 
data set was obtained in 1996/97 at 820 GeV and published in [1]. The results were obtained 
from two different samples, called A and B, covering higher > 12GeV^ and lower < 
12GeV^ regions, respectively. A comparison of the cross section measurement for sample A 
with the present analysis revealed a significant deviation, which required a dedicated study as 
described below. 
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5.1 Correction of the 820 GeV data 



For the analysis of the 820 GeV data an older DIS event simulation program version was used. 
If that program version is used without event weighting, the MC cross section agrees to better 
than 0.3% with the HECTOR calculation. However, in order to improve the statistical accuracy 
of the DIS MC event sample, the simulated events in [1] were generated with dependent 
weights for < 50 GeV^. This introduced a bias which depends only on Q^. The cross section 
as obtained in [1] can be corrected by a factor defined as 

^weighted , 

Here Nmc(Q^) and Nl^^^""^{Q^) are the sums of weights of events generated for a given bin 
in unweighted and weighted mode, respectively. The dependence of the correction factor on 
is empirically parameterised by 



^2^ 



Co fore2>50GeV2 
c(.Q ) = \ , , , , (3) 

1 Co + ci • logio(2V50 GeV^) for < 50 GeV^ 



with Co = 1.027 and ci = 0.0352. This correction procedure introduces an additional uncorre- 
cted cross section uncertainty of 0.5% on the corrected published data. 

For a further cross check the 1996/97 data are reanalysed within the framework of the 
present analysis. The luminosity determination is improved, resulting in an additional nor- 
malisation shift by -1-0.5%, corresponding to one third of a standard deviation of the quoted 
normalisation uncertainty. No further significant deviation is observed between the published 
results and the reanalysis. 

Tables 8 and 9 present the 1996/97 data from [1] corrected for the weighting bias and 
the small luminosity shift. These tables therefore replace the previous data for > 12GeV^, 
sample A. The measurements of F2 for y < 0.6 are extracted from the corrected cross sections 
accounting for the Fi influence, as described in section 5.3. The sample 5 in [1], extending to 
lower values of < 12GeV^, was not aff"ected by the MC weighting problem and has been 
combined with further HI data as described in [2]. 



5.2 Combined cross sections 

The corrected 820 GeV data and the present 920 GeV data are shown in figure 4 and combined 
to determine a new HI measurement at 12GeV^ < < 150GeV^. The combination of the 
data sets is based on the prescription introduced in [27] and developed further in [2]. The combi- 
nation uses the reduced cross section data and consistently treats the correlated and uncorrected 
uncertainty information, which is given in tables 6-9. The data are averaged for the region of 
large x, defined as in [2] hy y < 0.35, but kept separate at higher y where the reduced cross 
section depends significantly on the centre of mass energy via the Fl term. Small differences 
in the measurement grid are corrected for by adjusting cross section values for points close in x 
such that they correspond to a common x value, prior to the averaging. 
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zuuu 


III ^ aCdlC 


79 


W. JW 


E'^ linearity 






Polar angle 6(> 


— U.4D 


r> no 


LAr hadronic scale 


— U.OD 




LAI noisc 


— u.zz 


U.Wt 


SpaCal hadronic scale 






yp background 


0.11 - 


-0.11 


Luminosity 


0.64 - 


-0.46 



Table 3: Shifts in the central values of the systematic uncertainties determined for the combi- 
nation of the 1996/97 data at Ep = 820 GeV and the 2000 data at Ep = 920 GeV. The value for 
each shift is given in units of the original uncertainty. Note that for the two data sets some of 
the correlated systematic error sources were defined differently. 

In the combination of the two data sets, assumptions have to be made about the relationship 
between systematic uncertainty sources, which may not be fully independent between the dif- 
ferent analyses of the 1996/97 and 2000 data. Reasons for correlations between data sets are 
the similarity in the calibration procedure and the detector setup, whereas uncorrelated effects 
include variations with time, beam conditions and changes in the analysis procedure. To esti- 
mate the sensitivity of the averaged result to the correlation assumptions, different averages of 
the 1996/97 and 2000 data are compared, obtained considering all possible assumptions as to 
which systematic error sources are taken to be correlated or uncorrelated between the data sets. 
The observed variation of the averaged cross sections is typically at the per-mille level, which 
is negligible compared with the total uncertainty. The variations of the total uncertainties never 
exceed 10% of the final uncertainty. The assumption that the systematic uncertainties on the 
two data sets are uncorrelated yields the largest uncertainty. Thus no correlations between the 
error sources of the 1996/97 and 2000 data are considered in the combination procedure. 

The 1996/97 and 2000 data sets are fully consistent, as determined in the averaging proce- 
dure, with^^Q(/ Hdof = 51.6/61. As a result of the averaging, the central values of the correlated 
systematic uncertainties are modified. The shifts of the central values in units of the original 
uncertainties are given in table 3. It is observed that none of the absolute values of the shifts 
exceeds one standard deviation. 

The combined measurement of the reduced ep scattering cross section and its uncertainties 
are listed in tables 10 and 11. The correlated uncertainties are given as 14 new sources, 5,, after 
diagonalisation of the error matrix as explained in [2]. The combined HI data on the inclusive 
cross section measurement have total uncertainties of 1.3 - 2% in most of the phase space. 



5.3 F2 and its derivatives 



The cross section data are used to extract F2 for y < 0.6, correcting for the small influence of 
Fl using 

(Tr 

Fl = . (4) 

1 - mih 
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The values used for R = Fl/(F2 - Fi) are quoted in the cross section tables and are taken 
from the NLO QCD fit introduced below. Figure 5 shows the structure function F2 at fixed 
as a function of x together with HI data from lower [2] and higher [9]. The data are well 
described by the NLO QCD fit. The rise of F2 towards low x is thus established at a much 
improved accuracy compared with the first observations [28,29]. There is no indication for a 
saturation of this behaviour in the Q^, x region of study. 

Figure 6 shows the measurement of the structure function F2 at fixed values of x as a func- 
tion of 2^, compared with the QCD fit described below. The strong rise with at low x is a 
consequence of the large gluon density in this region. The data are well described by the QCD 
fit for 2^ > 3.5 GeV^. Also the extrapolation of the fit down to 2^ = 1.5 GeV^ gives a reason- 
able description of the data. At the largest x value covered by the data presented here, - 0.1, 
the structure function becomes almost independent of as the result of a compensation of 
quark and gluon contributions to the In derivative of F2. 

The DGLAP evolution equations [30-34] determine the derivative (5^2/5 In 2^)j taken at 
fixed jc, with the dominant low x contribution arising from gluon splitting into a quark-anti- 
quark pair. The measurement of this derivative is a powerful constraint on the gluon distribution 
xg and the strong coupling constant [35]. A study of (dF2/d\n Qp-)^ at low x has been pre- 
sented previously by the HI Collaboration [1]. The method described there is used to determine 
this derivative for the new, combined F2 data set, including the low data [2]. The results are 
shown in figure 7 for different x values as a function of . The dependence of the derivative 
on 2^ is well reproduced by the QCD fit. 

The rise of the structure function F2{x, Q^) towards low x may be quantified by the derivative 
A = -{d In In x)q2 which is shown in figure 8. The result is more accurate than the previous 
measurement [6] and extends to lower Q^. Within the uncertainty of the data, the derivative is 
constant at small x < 0.01, i.e. F2 for fixed is consistent with a power law F2 oc x~'^. Small 
departures from this behaviour, as are inherent in the QCD fit, cannot be excluded either. The 
value of A increases from about 0.1 to 0.3 in the region covered, from about 1 to 100 GeV^. 
Data from [2] allow this measurement to be extended to values below the region of validity 
of the DGLAP evolution. 



6 QCD Analysis 

The neutral current cross section measurements presented here, together with the measurements 
at lower [2] and the previously published NC and CC data at higher [7-9], provide an 
accurate HI data set for the determination of the parton density functions (PDFs) of the proton. 
A new QCD analysis, referred to as HI PDF 2009, is performed, which supersedes the previous 
HI PDF 2000 fit [9], as it relies on the more accurate new data. It also uses a general variable 
flavour number scheme (VFNS) treatment [36] of the heavy quarks, unlike the former fit, which 
used a zero mass variable flavour number scheme ("massless" scheme). 
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6.1 Framework and settings 

The QCD analysis decomposes the structure functions into a set of parton densities: the gluon 
xg, the valence quarks, xm„ and xd^, and the combined anti-up type and anti-down type quarks, 
xU = xu + xc and xD = xd + xs + xb. These are parameterised^ at a starting scale Ql and 
are evolved using the DGLAP evolution equations. An adjustment of the parton distribution 
parameters is performed to best fit the measured cross sections. 

The analysis is performed at NLO within the MS renormalisation scheme. The program 
QCDNUM [37] is used to solve the evolution equations. From the evolved parton distributions, 
the structure functions are calculated in the VFNS scheme [36,38], using 0(a^J coefficient 
functions for the calculation of Fl. The factorisation and renormalisation scales are both set 
to Q^. The^^ function is defined as for the cross section averaging discussed in section 5.2 
and [2]. It is minimised using the MINUIT package. The correlations between data points 
caused by systematic uncertainties are taken into account, following the numerical method pre- 
sented in [39,40]. The correlations between the systematic error sources of the different high 
data sets are treated as described in table 2 in [9]. The error sources of the data presented here 
and of the lower data [2] are all taken to be uncorrected with any other source. The only 
exception is a 0.5% luminosity uncertainty, which is associated with the theoretical calculation 
of the Bethe-Heitler cross section and which is common to all HI cross section measurements. 

Following [41], the masses of the charm and beauty quarks are set to m,. = 1 .4 GeV and nih = 
4.75 GeV. The value of the strong coupling constant is taken to be a^iM-^) = 0.1 176 [42]. The 
starting scale Ql is chosen to be slightly below the charm threshold, Ql = 1-9 GeV^. Hence at 
the starting scale the anti-quark densities simplify to xU(x) = xu{x) and xD{x) = x{d(x) + s(x)). 
The anti-strange quark density, xs{x) = fsxD(x), is taken to be a constant fraction, = 0.31, 
of xD at the starting scale [8]. A cut Q^ > Q^^^^ = 3.5 GeV^ is applied in order to ensure that 
the data used in the fit correspond to a kinematic domain where leading twist perturbative QCD 
can be used to predict the cross sections. Variations around these central values are taken into 
account as model uncertainties, as described below. 

6.2 Parameterisation 

The parton distributions xP are parameterised at Ql using the general form 

xPix) = Apjc^^(l - xf' [\+Dpx + Epx^ + ...]. (5) 

The specific choice of the parameterisations is obtained as follows: the parameters D, E, . . . are 
added one-by-one in an iterative procedure and kept only, if they reduce the by more than 3 
units. Parameterisations obtained are discarded if the structure functions F2 and Fl are negative 
anywhere in the range 10"^ < jc < 1 for 2^ > Ql. Parameterisations leading to very low valence 
quark distributions as compared to the total sea-quark density xS{x) = 2x(U(x) + D(x)) at large 

^ The previous HI PDF 2000 fit used a very similar decomposition of the quark flavours, but instead of the 
valence quark distributions xm„ and xd,,, the combined up and down quark distributions xU and xD were used. 
These determine xuu = xU - xU and xdo = xD - xD assuming symmetry between the sea-quarks and anti-quarks 
for each flavour. 
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Data Set 


i^aia rum la 


^2 
/vunc 


low 




[2] 


JO 


Q 


medium Q 




this measurement 


QQ 

yy 


O i . J 


e^p NC high Q^, 94 


- 97 


[7] 


1 jU 


QO 6 


e*p CC high Q^, 94 


-97 


[V] 


ZJ 


Zl .z 


e NC high Q-^, 98 


- 99 


[8] 


1 'XCt 
Ijy 


110 1 
1 IZ. 1 




QQ 


LoJ 


Zo 


I 1 .J 


e^p NC high (2^ 99 


-00 


[9] 


147 


137.4 


e^p CC high 99 


-00 


[9] 


28 


31.1 



Table 4: For each data set used in the HI PDF 2009 fit, the number of data points is shown, 
along with the;^^^ contribution determined using the uncorrelated errors only ixinc)- 

X, which dramatically fail to describe vp and vd fixed target measurements of the structure 
function xFi, [43], are also not considered. Some parameterisations of the xm„ density involving 
a quadratic E or cubic F term in x in the expansion (5) are examples of these. They are, however, 
included in the parameterisation uncertainty discussed below. The same applies to fits where 
the parton distributions are negative at very high x, as happens for some parameterisations of 
the gluon density involving an £ or F term. Among the possible parameterisations the one with 
the lowest is selected. This procedure, at Q^ = 1.9 GeV^, leads to the following choice: 

xg{x) = AgX^'{\ - xf' [l + D^x] , (6) 

xulx) = A,,x^'-{\-xf"\ (7) 

xdXx) = Ad^x^'^H - xf'\ (8) 

xU{x) = AuX^o(i _ ^fo^ (9) 

xD{x) = A5jc^°(l -jcf". (10) 

The normalisation parameters and A^^ are not fitted, but are obtained from the other pa- 
rameters via the quark counting rules. Since the existing data have a limited sensitivity to the 
behaviour of the valence quark distributions at low x, it is assumed that = 5^/ . Similarly, 
the behaviour of the up and down anti-quarks at low x is assumed to be governed by the same 
power, Bfj = Bf). As in [9], the normalisations of the U and D distributions are related by 
Afj = A[,{1 - f,) which corresponds to the usual assumption that d/u 1 as x ^ 0. Finally, the 
normalisation Ag of the gluon distribution is derived from the momentum sum rule. The total 
number of free parameters is thus equal to ten. 



6.3 Fit results 

The central fit of this analysis as specified above has a^^ of 587 for 644 degrees of freedom. 
The value for each data set is given in table 4. No significant tension is observed between 
the fit results and the systematic uncertainties on the low, medium, and high data sets. 

The fitted parameters of the distributions at the starting scale are given in table 5. The 
resulting parton distributions, including the total sea-quark density, are shown at 2^ = 4 GeV^ 
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0.105* 


-0.177 


2A2 




xD 


0.152 


-0.177* 


3.42 





Table 5: Fitted parameters corresponding to the distributions xg{x), xu„{x), xd„{x), xU{x) and 
xD{x) at the starting scale 2^ = 1.9 GeV^ (see section 6.2). The symbol * indicates that the cor- 
responding parameter is not a free parameter of the fit, but is derived from the other parameters. 

in figure 9. The inner error band describes the experimental uncertainty, obtained from the 
criterion l^x^ = 1 using the Hessian method as described in [40] and the numerical algorithm 
presented in [44] . The middle error band represents the experimental and model uncertainties 
added in quadrature. The model uncertainties are obtained by varying: 

• the charm mass between 1 .38 GeV and 1 .47 GeV; 

• the bottom mass m^, between 4.3 GeV and 5.0 GeV; 

• the strange fraction from 0.25 to 0.40; 

• the value of 2^ . , from 2.25 GeV^ to 5.0 GeV^; 

• the starting scale down to 1.5 GeV^. 

The resulting model uncertainty at low x is dominated by the sensitivity of the fit to the 2q 
variation ^ . 

In addition to the model uncertainty a new error contribution is introduced resulting from 
the parameterisation choice. As explained in section 6.2, alternative parameterisations leading 
to good fit quality but peculiar behaviour at large x are used to estimate the parametrisation 
uncertainties. An envelope of these fit solutions is built which is added in quadrature to the 
contributions of the experimental and model uncertainties. 

The distributions of xu„{x), xd^ix), xg{x) and xS{x) are shown at the starting scale = 
1.9 GeV^ and at = lOGeV^ in figures 10 a)-b) and figures 10 c)-d), respectively, with both 
linear and logarithmic scales. A comparison of the PDFs at the starting scale with their be- 
haviour in the DIS region, here represented by = lOGeV^, illustrates the rather dramatic 
influence of the DGLAP evolution on the sea-quark and gluon densities. At = 1.9 GeV^ the 
sea-quark density rises towards low x, in contrast to the gluon distribution which has a valence 
quark-like shape. The evolution rapidly changes the low-jc behaviour of the gluon distri- 
bution, which starts to rise similarly to the sea-quark distribution towards low x. In contrast, 

^ An uncertainty contribution due to the go variation is calculated by allowing to be as low as 1.5 GeV^. 
This causes an increase of the x^, which is comparable to the x~ change in a massless fit, when is varied from 
4 GeV^ to 2 GeV^, as has conventionally been done in the past. Since the VFNS scheme implementation requires 
2o ^ "^c' ^ default of = 1.9 GeV^ is chosen and the calculated uncertainty is symmetrised. 
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the non-singlet valence quark distributions evolve very slowly, as expected. As is shown in 
figures 10 c)-d), xg is the dominating parton distribution at low x and higher Q^, here displayed 
for = lOGeVl 



7 Summary 

A new measurement is presented of the inclusive double differential cross section for deep 
inelastic positron-proton scattering, e'^p — > e'^X, in the region of small Bjorken x, 2 ■ 10~^ < x < 
0.1, and four- momentum transfer squared, 12GeV^ < < 150GeV^. The data, corresponding 
to an integrated luminosity of about 22pb~\ were obtained with the HI detector at the ep 
collider HERA at beam energies = 27.6 GeV and Ep = 920 GeV. A small bias in a similar 
previously published data set, taken at Ep = 820 GeV, is found and corrected. The two data sets 
are then combined and represent the most precise measurement in this kinematic region to date, 
with typical total uncertainties in the range of 1.3 - 2%. 

The kinematic range of the measurement corresponds to a wide range of inelasticity y, from 
0.005 to 0.6. The data are used to determine the structure function F2ix, Q^), which is observed 
to rise continuously towards low x at fixed Q^. 

The high precision of the data allows new measurements to be presented of the partial deriva- 
tives of F2{x, Q^) with respect to and x. For x < 0.01, the derivative of lnF2 with respect 
to \nx confirms a power law dependence of F2 oc x~'^ with A depending primarily on Q^. The 
derivative of F2 with respect to In is a measure for the gluon distribution xg at low x. 

An NLO QCD fit to the HI data alone, including the new measurements presented in this 
paper, is performed. The fit implements a variable flavour treatment of heavy quark threshold 
effects. This new HI PDF 2009 fit supersedes the HI PDF 2000 previously obtained and provides 
a new determination of the gluon and quark densities of the proton including experimental, 
model, and parameterisation uncertainties. 
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Table 6: Reduced cross section cr^ for each Q^, x bin from the 2000 data. F2 is extracted with 
R taken from the QCD fit described in this paper using (4). The uncertainties on cr^ are quoted 
in %. 5tot is the total uncertainty determined as the quadratic sum of all uncertainties, ^stat is 
the statistical uncertainty, ^uncor represents the uncorrelated systematic uncertainty, y^;, y^'^ , 
Jde^ jEhni, Tnoise^ T^h^ds' Jyp ^^c bin-to-bin correlated uncertainties on the cross section 
measurement due to uncertainties on the SpaCal electromagnetic energy scale and linearity, 
positron scattering angle, LAr calorimeter hadronic energy scale, LAr calorimeter noise, and 
the photoproduction background, respectively. The overall normalisation uncertainty of 1.2% 
is not included. 
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Table 7: Reduced cross section cr,- and Fj for each Q^, x bin from the 2000 data, continued from 
table 6. 
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Table 8: Corrected reduced cross section cr^ from HI 1996/97 data, presented in the same format 
as originally in [1], with uncertainties quoted in %. The correction procedure is described in 
the text. F2 is extracted for y < 0.6 with R given by the QCD fit described in this paper. The 
overall normalisation uncertainty of 1.7% is not included. 
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-0.70 


0.17 


1.34 


0.00 


25 


3.20 


10"-' 


0.086 


0.898 


0.235 


0.898 


2.78 


0.79 


1.54 


1.80 


-0.77 


-0.23 


0.92 


0.00 






10"-' 


























25 


8.00 


10"' 


0.034 


0.677 


0.241 


0.677 


2.52 


0.92 


1.62 


1.11 


-0.68 


-0.72 


0.84 


0.00 


25 


1.58 


10"- 


0.018 


0.559 


0.205 


0.559 


3.71 


0.85 


1.57 


1.36 


-0.88 


-2.44 


-1.42 


0.00 


25 


5.00 


10"- 


0.005 


0.457 


0.180 


0.457 


7.54 


1.28 


3.39 


0.99 


-0.68 


-3.28 


-5.62 


0.00 


35 


5.74 


10"" 


0.675 


1.512 


— 


— 


2.65 


1.36 


2.01 


-0.37 


-0.70 


0.64 


0.08 


-0.55 


35 


8.00 


10"" 


0.484 


1.390 


0.232 


1.440 


2.18 


0.88 


1.75 


-0.87 


-0.49 


0.19 


0.05 


-0.28 


35 


1.30 


10"' 


0.298 


1.212 


0.231 


1.226 


1.76 


0.77 


1.45 


-0.21 


-0.77 


0.01 


0.01 


-0.02 


35 


2.00 


10"' 


0.194 


1.058 


0.229 


1 .063 


1.76 


0.77 


1.47 


-0.55 


-0.55 


0.00 


0.00 


0.00 






10"' 


























35 


5.00 


10"' 


0.077 


0.843 


0.231 


0.843 


2.67 


0.94 


1.57 


1.47 


-0.69 


-0.18 


1.06 


0.00 


35 


8.00 


10"' 


0.048 


0.738 


0.228 


0.738 


2.40 


1.00 


1.61 


1.02 


-0.68 


-0.15 


0.82 


0.00 


35 


1.30 


10"- 


0.030 


0.642 


0.211 


0.642 


3.24 


1.15 


1.69 


1.31 


-0.68 


-1.84 


-0.87 


0.00 


35 


2.51 


10"- 


0.015 


0.538 


0.200 


0.538 


4.14 


1.14 


1.67 


1.46 


-0.88 


-2.80 


-1.50 


0.00 


35 


8.00 


10"- 


0.005 


0.424 


0.089 


0.424 


9.21 


1.84 


3.55 


0.69 


-0.81 


-3.16 


-7.59 


0.00 


45 


1.30 


10"' 


0.383 


1.321 


0.230 


1.348 


1.94 


0.92 


1.75 


-0.24 


-0.23 


0.11 


0.03 


-0.05 


45 


2.00 


10"' 


0.249 


1.141 


0.228 


1.149 


1.75 


0.88 


1.47 


-0.25 


-0.54 


0.00 


0.00 


0.00 


45 


3.20 


10 


0.156 


1.009 


0.225 


1 .01 1 


1 .81 


0.94 


1.48 


—0.17 


—0.66 


0.00 


0.00 


0.00 


45 


5.00 


10--^ 


0.099 


0.898 


0.226 


0.899 


2.81 


1.08 


1.58 


1.64 


-0.53 


0.26 


1.08 


0.00 


45 


8.00 


10--^ 


0.062 


0.766 


0.219 


0.766 


2.50 


1.15 


1.61 


1.24 


-0.54 


0.07 


0.72 


0.00 


45 


1.30 


10"- 


0.038 


0.669 


0.21 1 


0.669 


2.85 


1.28 


1.69 


1.49 


-0.81 


-0.85 


-0.51 


0.00 


45 


2.51 


10"- 


0.020 


0.541 


0.175 


0.541 


4.28 


1.25 


1.66 


1.36 


-0.54 


-2.99 


-1.70 


0.00 


45 


8.00 


10"- 


0.006 


0.408 


0.099 


0.408 


7.60 


2.06 


3.52 


1.27 


-0.86 


-2.98 


-5.47 


0.00 


60 


2.00 


10"-^ 


0.332 


1.285 


0.225 


1.303 


2.10 


1.03 


1.76 


0.29 


-0.62 


0.04 


0.03 


0.00 


60 


3.20 


lU 


0.208 


1.086 


0.220 


1 .091 


1 .93 


1.07 


1 .5 1 


—0.65 


—0.32 


0.00 


0.00 


0.00 


60 


5.00 


10"' 


0.133 


0.929 


0.216 


0.930 


3.08 


1.25 


1.65 


1.96 


-0.56 


0.34 


0.96 


0.00 


60 


8.00 


10"' 


0.083 


0.829 


0.210 


0.829 


2.77 


1.31 


1.69 


1.51 


-0.47 


0.24 


0.78 


0.00 


60 


1.30 


10"- 


0.051 


0.705 


0.199 


0.705 


2.94 


1.43 


1.74 


1.49 


-0.76 


-0.87 


-0.40 


0.00 


60 


2.00 


10"- 


0.033 


0.616 


0.184 


0.616 


4.07 


1.70 


1.88 


1.84 


-0.88 


-2.34 


-0.71 


0.00 


60 


3.98 


10"- 


0.017 


0.522 


0.132 


0.522 


4.71 


1.82 


1.95 


0.88 


-0.32 


-3.03 


-2.25 


0.00 


60 


1.30 


10"' 


0.005 


0.371 


0.057 


0.371 


9.42 


3.03 


3.88 


1.38 


-0.56 


-3.98 


-6.82 


0.00 


90 




10"' 


























90 


5.00 


10"' 


0.199 


1.031 


0.208 


1.035 


2.20 


1.28 


1.61 


-0.67 


-0.67 


0.00 


0.00 


0.00 


90 


8.00 


10"' 


0.124 


0.872 


0.202 


0.873 


3.35 


1.52 


1.76 


1.97 


-0.51 


0.83 


0.94 


0.00 


90 


1.30 


10"- 


0.076 


0.751 


0.190 


0.752 


2.75 


1.61 


1.84 


1.06 


-0.27 


-0.23 


0.58 


0.00 


90 


2.00 


10"- 


0.050 


0.638 


0.168 


0.638 


3.87 


1.83 


1.94 


1.74 


-0.82 


-2.02 


-0.33 


0.00 


90 


3.98 


10"- 


0.025 


0.522 


0.128 


0.522 


3.83 


1.97 


1.95 


1.23 


-0.53 


-2.28 


-0.08 


0.00 


90 


1.30 


10"! 


0.008 


0.350 


0.049 


0.350 


4.84 


3.35 


2.06 


0.89 


0.16 


-1.69 


-2.08 


0.00 


120 


5.00 


10"' 


0.266 


1.043 


0.201 


1.051 


3.69 


2.36 


2.20 


0.39 


-1.80 


0.00 


0.00 


0.00 


120 


8.00 


10"' 


0.166 


0.866 


0.192 


0.868 


3.00 


1.94 


1.90 


-0.41 


-1.31 


0.00 


0.00 


0.00 


120 


1.30 


10"- 


0.102 


0.768 


0.179 


0.768 


4.61 


2.10 


2.09 


2.77 


0.49 


1.44 


1.68 


0.00 


120 


2.00 


10"- 


0.066 


0.623 


0.164 


0.624 


3.95 


2.26 


2.15 


2.33 


-0.36 


-0.64 


-0.43 


0.00 


120 


3.20 


10"- 


0.041 


0.576 


0.139 


0.576 


5.80 


2.77 


2.51 


2.65 


-0.49 


-3.25 


-1.34 


0.00 


120 


6.31 


10"- 


0.021 


0.477 


0.090 


0.477 


4.95 


3.05 


2.71 


1.57 


-0.41 


-2.25 


-0.50 


0.00 


120 


2.00 


10"! 


0.007 


0.322 


0.034 


0.322 


10.59 


4.76 


4.82 


1.22 


0.37 


-3.21 


-7.37 


0.00 


150 


2.00 


10"- 


0.083 


0.732 


0.155 


0.732 


8.85 


4.45 


3.86 


5.11 


0.85 


2.89 


2.95 


0.00 


150 


3.20 


10"- 


0.052 


0.568 


0.129 


0.568 


7.99 


5.05 


4.10 


4.56 


-0.53 


-0.66 


0.68 


0.00 


150 


6.31 


10"- 


0.026 


0.431 


0.088 


0.431 


8.97 


5.68 


4.44 


4.69 


-1.27 


-2.19 


-0.52 


0.00 


150 


2.00 


10"! 


0.008 


0.305 


0.038 


0.305 


12.86 


8.39 


7.12 


3.05 


-0.72 


-3.74 


-4.50 


0.00 



Table 9: Corrected reduced cross section (Jr and F2, for ij < 0.6, from HI 1996/97 data, contin- 
ued from table 8. 
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V 


IJ 


tTj. 


R 


'^2 


(5jivc lot 




■Siive.uncor 


' 1 


'1 


ave 
'3 


'4 


ave 
' 5 


^avc 
'6 


' 1 




^avc 
'9 


' 10 


' 1 1 


' 12 


Vive 
' 13 


yive 


12 


1.61 


10-* 


0.826 


1.236 


— 


— 


5.72 


4.11 


3.78 


0.97 


0.00 


-0.06 


-0.15 


-0.04 


-0.02 


0.02 


-0.14 


0.23 


-0.47 


0.33 


-0.46 


-0.08 


-0.06 


12 


1.97 


10-* 


0.675 


1.278 


— 


— 


3.48 


0.88 


2.15 


1.02 


1.82 


1.37 


-0.17 


0.16 


-0.02 


0.07 


0.07 


0.22 


-0.40 


0.22 


-0.38 


-0.10 


-0.07 


12 


2.00 


lo-* 


0.591 


1.318 


0.202 


1.387 


2.92 


1.21 


0.71 


0.92 


1.44 


-1.80 


-0.21 


-0.48 


0.14 


0.15 


-0.19 


-0.03 


-0.10 


-0.09 


0.11 


0.04 


0.08 


12 


3.20 


lO-* 


0.369 


1.238 


0.205 


1.259 


1.57 


0.86 


0.68 


0.96 


0.24 


-0.31 


-0.19 


-0.23 


-0.06 


0.07 


-0.15 


0.03 


-0.15 


-0.05 


0.12 


0.04 


0.08 


12 


3.20 


10-* 


0.415 


1.234 


0.205 


1.262 


2.19 


0.57 


1.73 


0.98 


0.39 


0.26 


-0.20 


0.01 


-0.09 


0.00 


0.17 


-0.03 


-0.29 


0.19 


-0.30 


-0.04 


-0.10 


12 


5.00 


10-* 


0.236 


1.164 


0.209 


1.171 


1.40 


0.65 


0.67 


0.96 


0.03 


-0.03 


-0.22 


-0.10 


-0.10 


0.00 


0.02 


0.01 


-0.29 


0.05 


0.01 


0.03 


0.03 


12 


8.00 


10-* 


0.148 


1.064 


0.213 


1.067 


1.62 


0.94 


0.69 


0.96 


0.00 


-0.03 


-0.30 


0.01 


-0.13 


-0.06 


-0.09 


0.10 


-0.41 


0.00 


0.19 


0.07 


0.09 


15 


2.01 


10-* 


0.827 


1.269 


_ 


_ 


4.98 


3.20 


3.61 


0.97 


0.00 


-0.06 


-0.15 


-0.05 


-0.03 


0.01 


-0.11 


0.20 


-0.47 


0.34 


-0.47 


-0.07 


-0.06 


15 


2.46 


10-* 


0.675 


1.369 


— 


— 


2.94 


0.92 


2.17 


0.99 


0.93 


0.67 


-0.18 


-0.02 


-0.04 


0.01 


0.00 


0.13 


-0.52 


0.42 


-0.54 


-0.09 


-0.03 


15 


3.20 


10-* 


0.462 


1.329 


0.214 


1.369 


1.91 


0.84 


0.70 


0.94 


0.70 


-0.87 


-0.06 


-0.44 


0.07 


0.11 


-0.15 


0.00 


-0.16 


-0.07 


0.11 


0.04 


0.08 


15 


3.20 


10-* 


0.519 


1.307 


0.214 


1.359 


2.51 


0.67 


1.97 


0.98 


0.65 


0.46 


-0.15 


0.03 


-0.05 


-0.03 


-0.15 


0.17 


-0.33 


0.20 


-0.32 


-0.04 


-0.11 


15 


5.00 


10-* 


0.295 


1.225 


0.217 


1.238 


1.34 


0.55 


0.66 


0.96 


0.07 


-0.07 


-0.16 


-0.14 


-0.08 


0.04 


-0.12 


0.04 


-0.20 


0.00 


0.04 


0.03 


0.04 


15 


8.00 


10-* 


0.185 


1.119 


0.220 


1.123 


1.35 


0.55 


0.64 


0.96 


0.01 


-0.03 


-0.23 


-0.02 


-0.10 


-0.01 


-0.06 


0.05 


-0.30 


0.04 


0.04 


0.04 


0.03 


15 


1.30 


10--' 


0.114 


0.992 


0.224 


0.993 


1.43 


0.59 


0.67 


0.96 


-0.01 


-0.03 


-0.29 


0.04 


-0.14 


-0.13 


0.04 


-0.01 


-0.39 


-0.07 


0.21 


-0.05 


0.00 


15 


2.00 


10--' 


0.074 


0.883 


0.228 


0.883 


1.49 


0.61 


0.67 


0.97 


0.00 


-0.03 


-0.24 


0.06 


-0.39 


-0.34 


-0.24 


-0.07 


0.10 


-0.27 


0.01 


-0.05 


0.00 


15 


3.20 


10--' 


0.046 


0.797 


0.233 


0.797 


1.46 


0.66 


0.69 


0.97 


0.00 


-0.02 


-0.25 


-0.01 


-0.22 


-0.21 


-0.27 


-0.10 


0.09 


-0.24 


0.02 


-0.05 


-0.01 


15 


5.00 


10--' 


0.030 


0.728 


0.237 


0.728 


1.49 


0.70 


0.70 


0.97 


0.00 


-0.02 


-0.23 


-0.04 


-0.28 


-0.26 


-0.21 


-0.06 


0.14 


-0.20 


0.03 


-0.01 


-0.01 


15 


1.00 


10-- 


0.015 


0.609 


0.242 


0.609 


1.38 


0.59 


0.65 


0.96 


0.00 


-0.02 


-0.20 


0.00 


-0.12 


-0.17 


-0.21 


-0.10 


0.17 


-0.13 


0.04 


0.08 


-0.01 


15 


2.51 


10-- 


0.006 


0.505 


0.224 


0.505 


2.63 


0.72 


1.79 


0.96 


-0.01 


-0.03 


0.20 


0.29 


-0.12 


0.98 


-0.02 


0.43 


0.61 


0.67 


0.26 


0.26 


-0.08 


20 


2.68 


10-* 


0.826 


1.334 


— 


— 


5.05 


3.25 


3.67 


0.97 


0.00 


-0.06 


-0.15 


-0.03 


-0.02 


0.03 


-0.12 


0.23 


-0.46 


0.32 


-0.46 


-0.08 


-0.06 


20 


3.28 


10-* 


0.675 


1.408 


— 


— 


2.79 


1.04 


2.10 


0.99 


0.78 


0.56 


-0.17 


0.10 


-0.03 


0.08 


0.07 


0.20 


-0.36 


0.18 


-0.33 


-0.09 


-0.09 


20 


5.00 


10-* 


0.394 


1.319 


0.225 


1.348 


1.50 


0.74 


0.70 


0.96 


0.26 


-0.33 


-0.12 


-0.20 


-0.01 


0.06 


-0.16 


0.00 


-0.12 


-0.05 


0.11 


0.04 


0.08 


20 


5.00 


10-* 


0.443 


1.321 


0.225 


1.359 


2.16 


0.64 


1.71 


0.97 


0.26 


0.16 


-0.13 


0.04 


-0.05 


-0.03 


-0.22 


0.21 


-0.27 


0.10 


-0.23 


-0.04 


-0.14 


20 


8.00 


10-* 


0.246 


1.209 


0.227 


1.218 


1.30 


0.49 


0.64 


0.96 


0.02 


-0.05 


-0.15 


-0.05 


-0.07 


-0.02 


-0.13 


0.06 


-0.23 


0.04 


0.02 


0.03 


0.03 


20 


1.30 


10-' 


0.151 


1.087 


0.229 


1.089 


1.34 


0.51 


0.64 


0.96 


0.00 


-0.03 


-0.17 


-0.01 


-0.09 


-0.06 


-0.08 


0.08 


-0.34 


0.06 


0.05 


0.04 


0.03 


20 


2.00 


10-' 


0.098 


0.979 


0.231 


0.980 


1.44 


0.55 


0.66 


0.96 


-0.01 


-0.03 


-0.22 


0.03 


-0.13 


-0.16 


-0.04 


0.06 


-0.49 


-0.04 


0.23 


-0.02 


0.00 


20 


3.20 


10-' 


0.062 


0.863 


0.233 


0.864 


1.38 


0.57 


0.66 


0.96 


0.00 


-0.02 


-0.12 


0.00 


-0.23 


-0.23 


-0.17 


-0.06 


0.12 


-0.21 


0.02 


-0.02 


0.01 


20 


5.00 


10-' 


0.039 


0.774 


0.235 


0.774 


1.39 


0.61 


0.68 


0.96 


0.00 


-0.02 


-0.22 


-0.01 


-0.11 


-0.22 


-0.14 


-0.04 


0.13 


-0.16 


0.02 


-0.01 


0.01 


20 


1.00 


10-2 


0.020 


0.653 


0.235 


0.653 


1.37 


0.49 


0.63 


0.96 


-0.01 


-0.01 


-0.19 


-0.03 


-0.23 


-0.17 


-0.21 


-0.20 


0.12 


-0.28 


-0.03 


0.06 


0.02 


20 


2.51 


10-2 


0.008 


0.524 


0.207 


0.524 


2.37 


0.56 


1.78 


0.96 


-0.01 


-0.02 


0.00 


0.23 


-0.39 


0.82 


-0.16 


0.04 


0.40 


0.08 


0.06 


0.31 


-0.04 


25 


3.35 


10-* 


0.827 


1.407 


_ 


_ 


5.80 


4.09 


3.92 


0.97 


0.00 


-0.06 


-0.15 


-0.04 


-0.03 


0.02 


-0.12 


0.20 


-0.45 


0.32 


-0.45 


-0.07 


-0.06 


25 


4.10 


10-* 


0.675 


1.408 








2.79 


1.15 


2.10 


0.99 


0.72 


0.51 


-0.14 


0.18 


0.00 


0.14 


0.03 


0.33 


-0.28 


0.03 


-0.22 


-0.10 


-0.12 


25 


5.00 


10-* 


0.492 


1.405 


0.229 


1.457 


1.79 


0.94 


0.75 


0.95 


0.52 


-0.65 


-0.12 


-0.25 


0.11 


0.11 


-0.19 


-0.05 


-0.08 


-0.08 


0.09 


0.03 


0.08 


25 


5.00 


10-* 


0.553 


1.376 


0.229 


1.445 


2.46 


1.04 


1.88 


0.97 


0.32 


0.20 


-0.18 


0.00 


-0.06 


0.00 


0.04 


0.08 


-0.35 


0.23 


-0.35 


-0.06 


-0.09 


25 


8.00 


10-* 


0.308 


1.279 


0.230 


1.294 


1.33 


0.54 


0.67 


0.96 


0.04 


-0.05 


-0.13 


-0.09 
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Table 10: Combined reduced cross section measurement from HI. The uncertainties are quoted in %. ^ave.tot is the total uncertainty deter- 
mined as the quadratic sum of all uncertainties, ^ave.stat is the statistical uncertainty, ^ave.uncor represents the uncorrelated systematic uncer- 
tainty, y^'^^ are the bin-to-bin correlated systematic uncertainties on the cross section measurement after diagonalisation. All normalisation 
uncertainties are included, except for the correlated 0.5% common to all HI cross section measurements. 
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041 
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0.45 
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0.24 


1.11 


-0.40 



Table 11: Combined cross section measurement from HI, continued from table 10. 




Figure 1: Distribution of events requiring < 0.6 and «/s > 0.005: the polar angle a) and the 
energy b) of the scattered positron, the transverse momentum ratio P^/P^ (data outside the anal- 
ysis selection indicated by the vertical line shown as open symbols) c) and E - d). The curves 
represent the MC simulation normalised to the luminosity, including a very small contribution 
from photoproduction. The narrow bands illustrate the correlated systematic uncertainties of 
the measurement. 
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Figure 2: Bjorken-x and distributions of events requiring ye < 0.6 and y^. > 0.005: Xe 
and Ql reconstructed using the electron method a)-b), and x-^ and reconstructed using the 
D method c)-d). The curves represent the MC simulation normalised to the luminosity. The 
narrow bands illustrate the correlated systematic uncertainties of the measurement. The small 
photoproduction background is drawn shaded. 
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Figure 3: Comparison of reduced cross sections as obtained from the Ep = 920 GeV data 
with the electron (closed circles) and 2 (open squares) reconstruction methods. The error bars 
represent the total measurement uncertainties. 
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Figure 4: Comparison of reduced cross sections as obtained from the Ep = 920 GeV (closed 
circles) and the corrected results from Ep = 820 GeV (open squares). The error bars represent 
the total measurement uncertainties. 
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Figure 5: Measurement of the structure function Fi at fixed as a function of x. The data of 
this measurement (closed circles) are complemented by the previously published data at low 
(open circles) [2] and high (open boxes) [9]. The error bars represent the total measurement 
uncertainties. The curve represents the QCD fit described in this paper. 
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Figure 6: Measurement of the structure function F2 as a function of at various values of x. 
The data of this measurement (closed circles) are complemented by the previously published 
data at low (open circles) [2] and high (open boxes) [9]. The error bars represent the total 
measurement uncertainties. The solid curve represents the QCD fit described in this paper for 
> 3.5 GeV^, which is also shown extrapolated down to = 1.5 GeV^ (dashed). 
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Figure 7: Logarithmic derivative of the structure function F2 as a function of at various 
values of x. The data of this measurement (closed circles) are complemented with the published 
data at lower (open circles) [2]. The error bars represent the total measurement uncertainties. 
The solid curve represents the prediction of the QCD fit for > 3.5 GeV^, which is also shown 
extrapolated down to = 1.5 GeV^ (dashed). 
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Figure 8: Measurement of the function A{x, Q^), defined as the negative logarithmic x derivative 
of lnF2 as a function of x at various values of Q^. The data of this measurement (closed 
circles) are complemented with the published data at lower (open circles) [2]. The error bars 
represent the total measurement uncertainties. The solid curve represents the prediction of the 
QCD fit for > 3.5 GeV^, which is also shown extrapolated down to = 1.5 GeV^ (dashed). 
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H1 PDF 2009, Q^ = 4GeV^ 




Figure 9: Parton distributions as determined by the HIPDF 2009 QCD fit at = 4GeV^. 
Shown are the combined up and down quark distributions, xU = x{u + c) and xD = x{d + s), 
their anti-quark counter parts, xU and xD, the valence quark distributions, xuv and xdv, the total 
sea distribution, xS = 2x{U + D), and the gluon distribution, xg. The inner error bands show 
the experimental uncertainty, the middle error bands include the theoretical model uncertainties 
of the fit assumptions, and the outer error band represents the total uncertainty including the 
parameterisation uncertainty. 
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Figure 10: Parton distributions as determined by the HI PDF 2009 QCD fit at = 1.9 GeV^ 
a)-b) and at = 10 GeV^ c)-d). In a) and c) (linear vertical scale), the gluon and sea-quark den- 
sities are downscaled by a factor 0.05. The inner error bands show the experimental uncertainty, 
the middle error bands include the theoretical model uncertainties of the fit assumptions, and 
the outer error band represents the total uncertainty including the parameterisation uncertainty, 
see text. 
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